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Abstract 

We propose a relativistically covariant model of interacting dark energy based on the 
principle of least action. The cosmological term A in the gravitational Lagrangian 
is a function of the trace of the energy-momentum tensor T. We find that the A(T) 
gravity is more general than the Palatini f(R) gravity, and reduces to the latter if 
we neglect the pressure of matter. We show that recent cosmological data favor a 
variable cosmological constant and are consistent with the A(T) gravity, without 
knowing the specific function A(T). 
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1 Introduction 



The most accepted explanation of the observed cosmic acceleration [1] is that 
the universe is dominated by an exotic component with large negative pres- 
sure, called dark energy [2] . The simplest way of introducing dark energy into 
general relativity is to add a cosmological constant A to the Riemann curvature 
scalar R in the Lagrangian: 

S = S m -^- I d'x^lR + 2A], (1) 



which leads to Friedmann and Lemaitre cosmological models [3]. The A term is 
spatially uniform and time independent (explicitly), as required by the princi- 
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pie of general covariance, and can be viewed as a representation of a relativistic 
ideal fluid obeying the equation of state p = — e. Moreover, this fluid obeys 
the equation of continuity that does not depend on matter energy density. 
Therefore, such a form of dark energy is said to be non-interacting. 



An interaction between ordinary matter and dark energy can be introduced in 
the form of a time dependent cosmological constant A(t) [4]. A decreasing A 
gives a positive rate of the entropy production and can be responsible for the 
observed large entropy of the universe. The cosmological constant can depend 
on the cosmic time through the scale factor a or the Hubble parameter H, and 
dimensional arguments lead to A oc a" 2 or A oc H 2 , respectively [5,6]. This 
dependence can, of course, have a different form and involve other quantities 
such as temperature. Phenomenological models with a variable cosmological 
constant are listed and reviewed in [7]. However, they all lack a Lagrangian 
description, i.e. the interaction between matter and dark energy cannot be 
derived in these models from the principle of least action in a relativistically 
covariant form. 



To preserve the covariance, a variable cosmological constant must depend only 
on relativistic invariants. A possible choice is A as a function of R. Such mod- 
els are referred to as f(R) gravity and can explain current cosmic acceleration, 
e.g., for A(R) oc R~ x [8,9]. In these models, a Legendre-Helmholtz transfor- 
mation of the Lagrangian or a conformal transformation of the metric bring 
the gravitational field equations of f(R) gravity into the form of the Einstein 
equations of general relativity with an additional scalar field [10]. A similar 
equivalence can be proven for a Lagrangian depending on invariants composed 
of the Ricci tensor R^ v [11]. 



In this Letter, we propose another choice in which the cosmological constant 
is a function of the trace of the energy-momentum tensor, and dub such 
a model A(T) gravity. This choice has an advantage over f(R) and f{R lxv ) 
gravity theories because we do not need to transform (conformally) from one 
metric to another, and there is no question about which metric is physical 
since we use one and the same metric tensor [10]. In Sec. 2 we derive the field 
equations in A(T) gravity and apply them to the Robertson- Walker metric 
which describes a homogeneous and isotropic universe. In Sec. 3 we relate 
our model to f(R) gravity and some phenomenological models of interacting 
dark energy. In Sec. 4 we show that the cosmological data favor a variable 
cosmological constant. The results are summarized in Sec. 5. 
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2 The field equations and cosmology 



The A(T) gravity action is given by 

S = Sm '^cI d4x ^[R + 2A(T)], (2) 



and the dynamical energy-momentum tensor T pu is generated from the action 
for matter S m : 

5S m = ^-J'^x v ^T flu 5g^. (3) 

From the principle of least action, the variation of (2) with respect to the 
metric tensor g pv equals zero, 

5S m-^- c J d A x^(R, v - X -R 9ilv + 2K'{T)T, V - K{T)g^8g^ 
_1 f d A x^g-M{T)g^8T^ = 0, (4) 

KC J 

where the prime denotes the differentiation with respect to T. For the variation 
of T pv we write a general expression, 

ST, U = N^ScT = -N^5g pa , (5) 



where the tensor N puprj depends on the kind of matter described by T pu . 
Eqs. (4) and (5) yield the Einstein field equations, 1 

iV - \Rg,u = [« - 2A'(T)]T, U + A(T)^„ - 2A\T)N" piu/ . (6) 



For a massless scalar field, T pu = ^d^d^cj) — ^g pu d p (f)d p (f), we find 



N pvp a = ^g P ( M gv)*d x <f)d x <i) - ^g pv d p <pd a <p, (7) 
N p ppu = \g,vd p <pdo<P - 8,<pd u <P = -2T, U + ±Tg, v , (8) 



1 The cosmological term A(T) modifies the gravitational interaction between matter 
and curvature, replacing k by the running gravitational coupling parameter which 
is a function of the tensors T pv and N p pfiu . 



3 



where the brackets denote symmetrization. In this case, Eq. (6) becomes 

iV - \Rg^ = [« + 2A'(T)]T^ + [A(T) - TA'(T)]^. (9) 



For the electromagnetic field, T pv = \g llv F pa F pc ' — F lip F v p , we obtain 



N^vpa = F^F^ + -g^ v F pX F a x - -g p ^g u)a F KX F KX , (10) 

N p ppu — Fp P F v p — -gp U F p<J F p = —T pu . (11) 

In this case, Eq. (6) reads simply 

iV - = kT^ + A(T)g pu . (12) 



The term A(T) does not affect the electromagnetic field since T = 0. For a 
perfect fluid, T pv = (e + £>)?v ~~ w e derive 

A^i/p* = P9p(n9u)a ~ (e + p)(g P ( M u u )U a + g a (nU u )U p ), (13) 

= - 2(e + = —1Tp V - pgp U . (14) 

Unlike for the previous cases, the tensor iV** cannot be now expressed as a 
linear combination of Tp U and Tg^. 2 The field equations (6) give 

R; - ^ = [« + 2A'(T)](e + p)«X + [A(T) - M^. (15) 



The Bianchi identity requires the right-hand side of (6) to be covariantly 
conserved, which leads to 



= {2A"(T)T'%„ - A'(T)Tp + 2A"{T)T* 
+2A'(T)iV^^}[K-2A'(T)]- 1 . (16) 



2 Contracting Eq. (6) with the metric tensor gives —R = [k — 2A' (T)]T + 4A(T) — 
2A' \T)N P a . The impossibility of expressing the tensor N p p ^ u as a function of 
and Ttjyu, means that the scalar N p pa a is not a function of T only. For some special 
cases, such as a massless scalar field or a pressureless dust, the quantity N p pa a is a 
definite function of T and so is R. Therefore, if we know the function A(T), we obtain 
an algebraic equation for T = T(R) and the A(T) gravity reduces to the Palatini 
f(R) gravity formulated in the Einstein frame (see the end of Sec. 3) [10,12]. 
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Eq. (16) means that the energy-momentum tensor is not, in general, covari- 
antly conserved. We can write the field equations (6) as 



R»„--Rg^ = K(T™+Tt 1/ ), (17) 



where T™ = T^ u , which defines the dark energy-momentum tensor, 3 



1%, = -\A(T)g, v - 2A> {T){T, V + N> 



From Eq. (17) it follows that matter and dark energy form together a system 
that has a conserved four-momentum. Therefore, we can speak of an interac- 
tion between matter and dark energy that together form a closed system. For 
a perfect fluid, we can write Eq. (18) as 

T A = ( Ca + P^) u v u » - VhQuv, (19) 



from which we obtain the expressions for the dark energy density and 
pressure p^. 

2A'(T) , A(T) . , 

CA = —{tm + Pm) + (20) 

K K 

P , = (21) 

K 

where T = e — 3p. 

We now apply Eq. (15) to a flat [13], homogeneous and isotropic universe filled 
with a perfect fluid, described by the Robertson- Walker metric, 

ds 2 = c 2 dt 2 -a 2 (t)(dr 2 + r 2 dQ 2 ). (22) 
The first Friedmann equation gives the Hubble parameter, 



H= o_ = c (k + 2A')e m + 2A' Pm + A 



a 



Such defined dark energy-momentum tensor depends on the tensor T^ u . There- 
fore, the former represents both "pure" dark energy (vacuum) and the interaction 
between vacuum and matter. 
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where A = A(T) and A' = A'(T). The dot denotes the differentiation with 
respect to the cosmic time t. The second Friedmann equation reads 



a 2 + 2aa . /nA . 

C 2 a 2 = ~ K Pm + A, (24) 



which yields the deceleration parameter, 

= _aa = (k + 2A / )e m + (3k + 2A> m - 2A 
q a 2 2[(k + 2A')e m + 2A'p m + A] ' 1 j 



The density ratio parameter JIa = caAc, where e c = 3H 2 /(kc 2 ) is the critical 
energy density, is given by 4 

2A'(e m + Pm ) + A 

A (K + 2A')e m + 2A' Pm + A- 1 ) 



We restrict our attention to a universe filled with dust, for which p m = and 
T = e. Hereinafter, the prime denotes the differentiation with respect to the 
energy density e, and we write e instead of e m . Eqs. (23) and (25) become, 
respectively, 5 



H = c f + 2 f + A , (27) 
- (- + 2A ')-2A 



* 2[(/c + 2A')e + A] ' 
while Eq. (26) reduces to 

2A ' e + A 

n » = R«Ta' < 29) 



When the universe is dominated by matter, q ~ 1/2 and ~ 0, and when 
it is dominated by dark energy, g ~ — 1 and ~ 1. The deceleration-to- 
acceleration transition (q = 0) takes place when e = e t , where e t is determined 
by 

[k + 2A'(e t )]e t - 2A(e t ) = 0. (30) 



4 The €a defined in Eq. (20) satisfies e rn + eA = e c - Therefore, = £l m + f^A = 
(e m + €k)/e c = 1, as in general-relativistic cosmology for a flat universe. 

5 In the case of dust, Eq. (15) can be written as R^ v — \Rg ilu = Kj-T^ + A(T)g /Jtl/ , 
where the running gravitational coupling is given by n r = n + 2A'(T). 
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Applying the Bianchi identity to Eq. (17) leads to the equation of continuity, 
{e + e A y + 3H{e + e A + p A ) = 0, (31) 



which gives the time dependence of the matter energy density, 

e = -3He ± 2A> A . (32) 

k + 3A' + 2A"e 1 ; 

Using a = ao(l + z)~ l and integrating (32) gives the relation between the 
energy density and redshift, 

ln(1 + " ) = 3i7 . + 2A ' (33) 



where eo is the present energy density. Combining (27) with (33) yields the re- 
lation H = H(z) which is of observational interest. If A = const, we reproduce 
the AC DM relation H 2 {z) = y[«e (l + zf + A]. 

We can write the equation of continuity (31) as 



e + 3He = Q, (34) 
e A + 3H(e A + PA ) = -Q, (35) 

where Q is the rate of the energy transfer from dark energy to ordinary matter, 

A' + 2A"e , x 

Q = 3He — — — — . 36 

^ /t + 3A + 2A"e v ; 

The relative transfer rate T equals Q/e A , 

A' + 2A"e 

1 = 3KHe (« + 3A' + 2A"e)(2A'e + A)- (37) 

If A = then both Q and e A vanish and T is indeterminate. The non- 
dimensional interaction rate r = Y/H is given by 

3^(A + 2A"e) 

T (« + 3A' + 2A"e)(2A'e + A)' 1 ' 

The first and second law of thermodynamics for an isotropic and homogeneous 
universe can be written as (ea 3 )' + p(a 3 )' = TS, which with Eq. (34) yields 

TdS = a 3 Qdt. (39) 
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The quantity Q is related to the production of entropy in a growing universe, 
and thus must be positive. 



3 A(T) gravity and phenomenological models 



If A = const, we return to the AC DM cosmology. If A is proportional to 
T = e, A = riKe, where n ^ is a non-dimensional constant, we reproduce the 
model with A oc H 2 [6] : 

3nH 2 1 . , 

A= (TT^' r = rnw (40) 



where the Hubble parameter is given by 



H = M /,vr(rc+^). (11) 



We have a static case H — for n — — |, and no solution if n < — |. From 
non- negativity of Q and e it follows that n > 0. Eq. (32) and the definition of 
H yield the law for scaling of the energy density, 

eoca" 3 ^). (42) 



There is no deceleration-to-acceleration transition, however, since 

q = 2{l^n) = COnSt ' (43) 



To explain the observed transition, we need to break the proportionality A oc e. 
The simplest choice is A = nne + Ao, where Ao is a constant 6 . The scaling 
law (42) does not change, while the deceleration parameter becomes 

K€ ~ 2A ° (44) 



2[(l + 3n)«e + Ao]' 



The transition occurs at 

e t = 2 A (45) 

K 



This form can be viewed as the linear approximation of the real function A(e). 
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from which we obtain, using Eq. (33), the following relation: 



Zt = 



2An\ ( 3(l + 2n) ) 



€ Q K 



1. 



(46) 



By means of Eq. (29), we derive 



z t 



1 + 3n)Q^ — 3n 



1 — ^A,0 



( l+3n ) 
\ 3(l + 2n) / 



" 1, 



(47) 



where the subscript refers to the present time. The Hubble parameter as a 
function of the redshift can be found using Eqs. (27) and (42): 



H(z) = 



V3 



{ 3(l + 2n) \ I 1 / 2 

[1 + 3n)/te (l + z)v-i+srJ + A 



(48) 



By means of Eq. (20) and the definitions of Q m ,o and Q\,o, we bring Eq. (48) 
to the standard form, 



H(z) = H 



( 3(l + 2n) ^ 

(1 + 3n)Q mfi (l + z^^+^1 - 3nQ mfl + Q Afi 



1/2 



(49) 



For n = 0, Eqs. (47) and (49) reduce to their ACDM expressions. 
The relation A oc H 2 holds for a larger class of functions A(e) obeying 

(k + 2A')e + (1 - AT^A = 0, (50) 



where the constant of proportionality k is defined such that A 
The solution of (50) is given by 

kne i-fc 



3kH 2 /c 2 . 7 
(51) 



Hereinafter, C denotes a constant. The singular case A = H 2 /c 2 can be ob- 
tained from A = — |elne + Ce. The general condition for the constancy of r 
leads to a more complicated differential equation for A. 

Another phenomenological condition A oc a~ m , where m = const, can be 
found by combining (32) with the identity e = (de/dA)(dA/da)Ha. We find the 

7 The constants n and k are related by k = jrk^- 
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second-order differential equation for A(e) which is equivalent to this condition: 
mA(K + 3A' + 2A"e) - 3eA'(« + 2A') = 0. (52) 



This equation has a solution linear in e, A = 3 ( 3 m ™2) Ke ' exce Pt for the dimen- 
sionally favored case m = 2 [5] . 

The condition for no interaction between matter and dark energy follows 
from (32) or (36): 

A' + 2A"e = 0. (53) 



The solution of this equation is given by 

A = Ce 1/2 + A . (54) 



In this case, the deceleration parameter vanishes when 

Ke - Cyi - 2A = 0, (55) 



which has one physical solution for e. 

Finally, we examine the case in which A is a linear function of q. The relation 
A oc q gives a nonlinear differential equation for the function A(e). The Palatini 
f(R) gravity in the Einstein conformal frame predicts 

A = ^(l-2g), (56) 



which is linear in q but also includes the dimensional factor H 2 [14]. Using 
Eqs. (27) and (28) we can show that the same relation (56) is satisfied in the 
A(T) gravity. Therefore, we conclude that the A(T) gravity and the Palatini 
f(R) gravity are cosmologically equivalent, although this is true only for mat- 
ter without pressure. For more general cases, including matter described by 
both e and p, the A(T) gravity is broader than the Palatini f(R) gravity. In 
the latter, R is solely a function of T [8], which is generally not true for the 
A(T) gravity. 8 



See the second footnote in Sec. 2. 
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4 Observational constraints 



Recent WMAP cosmological data [15] indicate that matter in the universe 
scales almost like dust, ea 3 ~ const. If we define e as the deviation from the 
factor 3 in the scaling law, 

e ~ a-^~ £ \ (57) 



then the data require < e < 0.1 [16]. Assuming A = uke + A gives, using 
Eq. (42) and (57), the equivalent condition 

< n < 0.037. (58) 



The current value of the dark energy density parameter is Qa,o — 0.7li !o5 [13], 
which in the ACDM model gives Zt = 0.70^'^. This value of Zt slightly 
overlaps with the observed z t = 0.46 ± 0.13 [13]. For the largest n compatible 
with observations, n = 0.037, we use Eq. (47) to find 

z t = 0.70+°;}°. (59) 



This value slightly overlaps with the observed z t as well. Consequently, the 
A(e) = n/te + A model is consistent with astronomical observations, provided 
n <C 1, which means that the interaction between dark energy and ordinary 
matter is relatively weak [14]. 

We now proceed to show that the cosmological data favor a variable cosmo- 
logical constant. If we define 

a= — , P = — , 7 = , (60) 

K6q K K 



then we can express fi mj0 and q as 

a,„ = (i + a + 2, 9) -, * = ^ a a + + ll y (61) 



Reversing these equations yields 

^ = 1 - 2q p = 2 + 2q - 3^ m;0 2 
3f2 m n 6Jl m n 



In the ACDM model, (3 = and q = 3fi mj0 /2 - 1. 
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The values for the present matter density parameter and deceleration parame- 
ter from the SNIa gold sample data are Vl mfi = 0.29t°;o3 and q = -0.74±0.18, 
respectively [13]. They give 



a = 2.8512;??, ( 63 ) 
(3 = -0.20±g;l|. (64) 

The result of (64) favors a variable cosmological constant ((3 ^ 0), although 
does not exclude a possibility of A being a constant ((3 = 0). 9 The redshift of 
the deceleration-to-acceleration transition can be found from equations (30) 
and (33), and depends on a particular form of the function A(e). 

Another non-dimensional parameter measured by cosmologists is the current 
value for the redshift derivative of the deceleration parameter, (dq/dz) . We 
use the identity dq/dz = (dq/de)(dz/de)~ 1 , Eq. (33) written as 



dz 
de 



1 + 3/3 + 27 

36o(l + 2/3)' (65) 



and the derivative of Eq. (28) with respect to e, to obtain 



dq 

dz 



3(1 + 2(3) (2a -3(3 + 3a(3 - 6(3 2 + 4q 7 ) 

2(1 + 3(3 + 2^(1 + a + 2[3) 2 • ( j 



Let us assume (3 — 0. In this case, Eq. (66) becomes 
3a 



dq 

dz 



(1 + a) 2 ' 



(67) 



and is independent of 7. Substituting here a from (63) yields 



dq 

dz 



— u -°'-0.07' 



(68) 



which does not overlap with the observed 1.59 + 0.63 [13]. Therefore, the SNIa 
gold sample data for (dq/dz) support (3 7^ 0. We emphasize that the above 
derivation is general, since we did not use the specific function A(e). 



9 In the case A(e) = n/ce + Ao, we have (3 = n. Since n cannot be negative, we 
are left with the condition n < 0.06 which is on the order of the previously found 
inequality (58). 
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5 Summary 



The presented A(T) gravity is a relativistically covariant model that predicts 
a variable cosmo logical constant. It derives the field equations and the in- 
teraction between matter and dark energy from the principle of least action, 
generating this interaction from the cosmological term A which is a function 
of the trace of the energy-momentum tensor T. We showed that the A(T) 
gravity and the Palatini f(R) gravity give the same dynamics of the universe 
filled with dust, i.e. they are cosmologically equivalent if T = e. For more gen- 
eral cases, the A(T) gravity is broader. We also showed that cosmological data 
favor the cosmological constant which varies with time. Further astronomical 
observations should provide insight into the form of the function A(T), de- 
scribing how matter interacts with dark energy and how dark energy affects 
the gravitational coupling between matter and spacetime. 



Acknowledgment 



The author would like to thank Tomi Koivisto for valuable comments. 



References 



[1] A. G. Riess et al, Astron. J. 116, 1009 (1998); S. Perlmutter et al, Astrophys. 
J. 517, 565 (1999); N. W. Ralverson et al, Astrophys. J. 568, 38 (2002); C. B. 
Netterfield et al., Astrophys. J. 571, 604 (2002); C. L. Bennett et al, Astrophys. 
J. Suppl. Ser. 148, 1 (2003). 

[2] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003); T. 
Padmanabhan, Phys. Rept. 380, 235 (2003); E. J. Copeland, M. Sami, and 
S. Tsujikawa, hep-th/06030"57| (2006). 

[3] J. E. Felten and R. Isaacman, Rev. Mod. Phys. 58, 689 (1986). 

[4] M. Bronstein, Phys. Z. Sowjetunion 3, 73 (1933); M. S. Berman, Phys. Rev. D 
43, 1075 (1991). 

[5] M. Ozer and M. O. Taha, Phys. Lett. B 171, 363 (1986); W. Chen and Y. S. 
Wu, Phys. Rev. D 41, 695 (1990). 

[6] J. C. Carvalho, J. A. S. Lima, and I. Waga, Phys. Rev. D 46, 2404 (1992); J. 
A. S. Lima and J. C. Carvalho, Gen. Relativ. Gravit. 26, 909 (1994). 

[7] J. M. Overduin and F. I. Cooperstock, Phys. Rev. D 58, 043506 (1998). 

[8] D. N. Vollick, Phys. Rev. D 68, 063510 (2003). 



13 



[9] S. M. Carroll, V. Duvvuri, M. Trodden, and M. S. Turner, Phys. Rev. D 70, 
043528 (2004). 



[10] G. Magnano, gr-qc/9511027 (1995). 

[11] A. Jakubiec and J. Kijowski, Phys. Rev. D 37, 1406 (1988). 
[12] N. J. Poplawski, Class. Quantum Grav. 23, 2011 ibid 4819 (2006). 
[13] A. G. Riess et al, Astrophys. J. 607, 665 (2004). 
[14] N. J. Poplawski, |gr-qc/0607124"l (2006). 

[15] D. N. Spergel et al, Astrophys. J. Suppl. Ser. 148, 175 (2003). 
[16] P. Wang and X. H. Meng, Class. Quantum Grav. 22, 283 (2005). 



14 



